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Abstract 

We analyze the possibilities of pairing between two different fermion species 
in asymmetric matter at low density. While the direct interaction allows 
pairing only for very small asymmetries, the pairing mediated by polarization 
effects is always possible, with a pronounced maximum at finite asymmetry. 
We present analytical results up to second order in the low-density parameter 
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The recent experimental achievement of trapping fermionic alkali atoms at very low 
density and temperature fl[] demands theoretical estimates of the characteristic size of the 
pairing gaps A that might be observable in those systems. The canonical case of an attractive 
s-wave interaction leads (including polarization effects) to the well-known low-density result 
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where p = k F /2m is the chemical potential and a < the s-wave scattering length. However, 
under certain circumstances direct s-wave pairing is not possible: in the case of repulsion, 
a > 0, or, e.g., in spin-polarized one-component Fermi systems. 

A similar, particular system has recently been advocated for experimental study. It is a 
spin-polarized alkali gas composed of two different hyperfme levels of 6 Li ||. In this specific 
environment an (attractive) s-wave interaction exists only between the atoms at different 
levels, whereas atoms at like levels can only interact via p-waves. This report is dedicated 
to the quantitative study of the principal possibilities of pairing in such a system, and in 
particular of the dependence of pairing on the asymmetry of the system, which evidently is an 
important experimental parameter that influences strongly the magnitude of the observable 
gaps. 

More generally, we will assume a fermionic system composed of two distinct species 1 
and 2 (carrying definite spin orientations) of equal mass m and with densities p± and p 2 , 
or equivalently total density p = p\ + pi and asymmetry a — (p% — p2)/(pi + P2)- We also 
introduce the notation ki = kp, (i = 1,2) for the Fermi momenta ki = (67r 2 pj) 1 / 3 , and 
Pi = kf/2m for the two chemical potentials at low density. 

For the sake of presentation we will for the moment consider an idealized system without 
direct interaction between like particles 1-1 and 2-2, whereas 1 and 2 are interacting via a 
potential V with a s-wave scattering length a. We are only interested in the situation at 
very low density, ki\a\ <C 1, where the pairing properties are completely determined by the 
scattering length, or, equivalently, the low- momentum s-wave T-matrix To = Aira/m. We 
will now analyze the pairing gaps generated by this interaction. 

In the case of attraction, a < 0, and for very small asymmetries clearly by far the 
dominant process is the pairing generated by the direct s-wave interaction between different 
species, see Fig. |l|(a). The BCS theory generalized to asymmetric matter 0-0] yields the 
basic coupled equations for the determination of the gap function A fc , total density p, and 
density difference Sp, 
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with the Fermi function f{E) = [1 + exp((3E)] 1 and 
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E± = E k ±5fi , E k = yjel + A 2 , e k = e k - // , e fc = /c 2 /2m . (5) 

The chemical potentials of the species 1 and 2 are /ii j2 = /i ± 5/i and V is the bare po- 
tential acting between them. Throughout this report we will only determine gaps at zero 
temperature, where one has f(E) = 6{—E) and therefore 

1 - f{E k ) - f(E+) = 9{E k ) , (6) 
f{E£) - f(E+) = 9{-E k ) , (7) 

i.e., the unpaired particles are concentrated in the energy interval — Se, \x + Se], Se = 



5fi 2 — A 2 , which does not contribute to the pairing interaction. This situation is illustrated 
in Fig. [| that sketches the BCS momentum distributions of the species 1 and 2, according 
to Eqs. (|3|) and (ffl), for a positive asymmetry. The excess of particles 1 is located in the 
interval around fi, Pauli-blocking the gap equation. This leads to a rapid decrease of the 
resulting gap when increasing the size 2Se of the interval, i.e., the asymmetry. 

In the weak-coupling case, A Se fi, which is adequate in the low-density limit, the 
momentum distributions of the two species are very sharp and one obtains from Eqs. @ 
and (|): 

So 3 Se _ ,„ s 
a = « -- « 1 , 8 
P 2 fi 

i.e., the asymmetry a is directly proportional to Se. Analyzing also the gap equation in the 
context of a weak-coupling approximation, one obtains for the dependence of the gap on the 
parameter Se |5|Q. 

Sfi + Se = const. = A ^> A 2 = A 2 , - 2A Se , (9) 

where A is the gap in symmetric matter of the same total density. Consequently the gap 
as a function of asymmetry is: 



A V 3A 1 ; 

The gap vanishes at a max = 3A /4/i, which at low density is indeed an extremely small 
number, cf. Eq. ([I]). Therefore, for very small asymmetries already, pairing generated by 
the direct interaction between different species becomes impossible. 

For larger asymmetry only p-w&ve pairing between like species can take place, which 
in leading order in density is given by the polarization diagram shown in Fig. 0(b). We 
discuss in the following the gap of species 1 mediated by the polarization interaction due 
to species 2. Clearly the results are invariant interchanging 1 and 2. Quantitatively the 
relevant interaction kernel reads at low density 



with the static Lindhard function (pertaining to the species 2) 
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The factor 1/2 in Eq. (TfT| ) corrects for the fact that conventionally the Lindhard function 
contains a factor two for the spin orientations, which is not present in our case. It should be 
noted that it is the absence of exchange diagrams that renders this low-density polarization 
interaction attractive in contrast to the case of one species, where the polarization effects 
reduce the s-wave BCS gap by a factor (4e) 1//3 in the low-density limit 0-ffl, see Eq. ([I]). It 
is also remarkable that the polarization interaction is always attractive, depending only on 
the square of the scattering length a [|l(J . 

Projecting out the L = 1 partial-wave interaction, one obtains in particular |]TT| . |T2" | 
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This function is normalized in symmetric matter, g(l) = 1, and plotted in Fig. 0(a). The 
numerical factor (21n2 — l)/5 in Eq. ( ]i~4| ) is often replaced by the approximate value 1/13. 



Using now the general result for the (angle-averaged) L-wave pairing gap [10,13,14 
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and considering that for a pure polarization interaction I\, the leading order in density is 
T L = T L , one obtains 
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Taking into account the dependence of this expression on the two Fermi momenta k\ and k 2 , 
the final result for the variation of the pairing gap with asymmetry a for fixed total density 
p is therefore 
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A(p,0) 
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with 

h(a) = 1 
and the density parameter 
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U ^ = 13 & 2 ' <* = ( 3 - 2 <>) 1/3 . (20) 

The function h(a) is displayed in Fig. [|(b), where a = 1 corresponds to pure 1-matter. One 
notes a maximum at a m 0.478 with an expansion h(a) ~ 0.465 — 1.343(a — 0.478) 2 . This 
means that the gap at this asymmetry is enhanced by a factor e 0A6u ~ io°- 2u compared to 
the symmetric case [h(0) = 0]. Evidently, in the low-density limit u — ■> oo this represents 
an enormous amplification at finite asymmetry. Around this peak, the variation of the 
gap with asymmetry is well described by a Gaussian with width a = 1/(1. 64-*/^). As an 
illustration, Fig. [3](c) shows the ratio, Eq. (|T8|), for a value of the density parameter u = 100. 
A peak of the order of 10 20 is observed, that becomes rapidly more pronounced and narrower 
with decreasing density p (increasing u), although of course at the same time the absolute 
magnitude of the gap decreases strongly with decreasing density. 

Let us now briefly discuss the next-to-leading-order effects, namely additional contribu- 
tions of order (/cjCt) 3 in the denominator of the exponent of Eq. (|17|) . There are two principal 
sources of such effects, which are shown diagramatically in Fig. [|. The first one, Fig. ^(a), 



is the direct p-wave interaction [15] between like species that we have neglected before. 



Parametrizing the low-density p-wave T-matrix in the standard form Ti(k, k') ~ ATia\kk' /m, 



where a x is the p-wave scattering length, leads to a BCS gap [15 
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The second type of third-order contributions are polarization effects involving the s-wave 
scattering length a. In contrast to the case of a one-component system, where there are 
several relevant diagrams, in the present two-component system only one diagram exists. It 
is drawn in Fig. f|(b). Unfortunately it can only be computed numerically, which will not 
be attempted here. 

Finally, to fourth order, there is a large number of diagrams contributing to the in- 
teraction kernel. Apart from that, at this order it is also necessary to take into account 
retardation effects, i.e., the energy dependences of gap equation, interaction kernel, and self- 
energy need to be considered |l(|. All this can only be done numerically and was performed 



in Ref. |12[ for the case of a one-component system. 

In any case, the existence of higher-order corrections will not alter the main conclusions 
drawn so far, namely the presence of a strongly peaked Gaussian variation of the gap with 
asymmetry. They will, however, shift this peak to a different density-dependent location, 
and also modify the absolute size of the gap. Apart from that, the perturbative approach 
that is followed here, is clearly limited to the low-density range kp\a\ < 1. For larger density, 
different theoretical methods have to be used |T7|, which is still a difficult field of current 
investigation. 

In conclusion, we have studied the possibility of pairing in asymmetric fermion matter 
composed of two distinct species. In the low-density limit, the direct s-wave interaction 
between different species produces a gap ~ exp[7i/2k F a] only for very small asymmetries of 
the order of A //i, Eq. (|T0|), whereas for larger asymmetries the polarization- induced p-wave 
attraction between two like species produces a much smaller gap ~ exp[— 13(7r/2A;^a) 2 ], 
which extends however in principle over the whole range of asymmetry. In practice a sharp 
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maximum at a ~ 0.478 (pi/p2 ~ 2.83, fci/fo ~ 1.41) appears. Explicit expressions for 
the variation of these gaps with asymmetry were given. Higher-order corrections will only 
quantitatively modify these particular features of pairing in the low-density regime. Clearly 
the experimental observation of both types of pairing is supposedly difficult, in the first case 
due to the nearly perfect symmetry that is required, in the second one due to the extremely 
small size of the resulting gap. 

We acknowledge useful discussions with U. Lombardo, A. Ramos, and P. Schuck. This 
work was supported in part by the programs "Estancias de cientfficos y tecnologos extran- 
jeros en Espaha," SGR98-11 (Generalitat de Catalunya), and DGICYT (Spain) No. PB98- 
1247. J. M. P. wishes to acknowledge support from a postgraduate fellowship of the Fundacio 
Universitaria Agusti Pedro i Pons. 



6 



REFERENCES 



[1] B. De Marco and D. S. Jin, Science 285, 1703 (1999). 

[2] L. P. Gorkov and T. K. Melik-Barkhudarov, Sov. Phys. JETP 13, 1018 (1961). 
[3] H. Heiselberg, C. J. Pethick, H. Smith, and L. Viverit, Phys. Rev. Lett. 85, 2418 (2000). 
[4] H.-J. Schulze, A. Polls, A. Ramos, Phys. Rev. C63, 044310 (2001). 
[5] H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. Hulet, Phys. Rev. Lett. 76, 10 
(1996); 

M. Houbiers, R. Ferwerda, H. T. C. Stoof, W. I. McAlexander, C. A. Sackett, and R. 

G. Hulet, Phys. Rev. A56, 4864 (1997). 
[6] A. Sedrakian, T. Aim, and U. Lombardo, Phys. Rev. C55, 582 (1997); 

A. Sedrakian and U. Lombardo, Phys. Rev. Lett. 84, 602 (2000). 
[7] A. I. Akhiezer, A. A. Isayev, S. V. Peletminsky, and A. A. Yatsenko, Phys. Rev. C63, 

021304 (2001). 

[8] M. Yu. Kagan and A. V. Chubukov, JETP Lett. 47, 614 (1988). 

[9] M. A. Baranov, Yu. Kagan, and M. Yu. Kagan, JETP Lett. 64, 301 (1996). 
[10] D. Fay and A. Layzer, Phys. Rev. Lett. 20, 187 (1968). 
[11] M. Yu. Kagan and A. V. Chubukov, JETP Lett. 50, 517 (1989). 

[12] D. V. Efremov, M. S. Marenko, M. A. Baranov, and M. Yu. Kagan, Sov. Phys. JETP 
90, 861 (2000). 

[13] P. W. Anderson and P. Morel, Phys. Rev. 123, 1911 (1961). 

[14] W. Kohn and J. M. Luttinger, Phys. Rev. Lett. 15, 524 (1965). 

[15] L. You and M. Marinescu, Phys. Rev. A60, 2324 (1999). 

[16] A. S. Alexandrov and A. A. Golubov, Phys. Rev. B45, 4769 (1992). 

[17] S. Babu and G. E. Brown, Ann. Phys. 78, 1 (1973); 

T. L. Ainsworth and K. S. Bedell, Phys. Rev. B35, 8425 (1987). 



7 



FIGURES 



FIG. 1. The two possible lowest-order pairing interactions: (a) Direct s-wave interaction 
between different species, (b) Polarization-induced p-wave interaction between like species. Vo and 
To are the s-wave (L = 0) bare potential and T-matrix between species 1 and 2, respectively. 

FIG. 2. BCS momentum distributions of the species 1 and 2 in asymmetric superfluid matter. 

FIG. 3. (a,b) The functions g and h, appearing in Eqs. ( |T5| ) and (|19|), respectively, (c) The 
variation of the gap with asymmetry, Eq. (18), for a density parameter u = 100. 

FIG. 4. Third-order diagrams: (a) Direct p-wave interaction, (b) Polarization contribution. 
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